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TECHNICAL NOTE 3448 


THEORETICAL ANALYSIS OF INCOMPRESSIBLE FLOW THROUGH A RADIAL-INLET 
CENTRIFUGAL IMPELLER AT VARIOUS WEIGHT FLOWS 
I - SOLUTION BY A MATRIX METHOD AND COMPARISON 
WITH AN APPROXIMATE METHOD 

By Vasily D. Prian, James J. Kramer, and Chung-Haa Wu. 


SUMMARY 

A method for the solution of the incompressible, nonvlscous flow 
through a centrifugal impeller. Including the Inlet region, is presented. 
Several numerical solutions are obtained for four weight flows through 
an impeller at one operating speed. The results are presented in a 
series of figures showing streamlines and relative velocity contours. 

A comparison is made with the results obtained by using a rapid approxi- 
mate method of analysis. 


INTRODUCTION 

In order to provide the fundamental information about the internal 
flow necessary for the rational design of efficient centrifugal compres- 
sors, two-dimensional solutions of the potential flow throu^ centrifugal 
conpressors on both blade-to-blade and meridional surfaces have been ob- 
tained by means of relaxation methods in references 1 to 4. In addition, 
a three-dimensional potential-flow solution was obtained by similar means 
in reference 5. However, all these solutions are for iinpellers with in- 
ducer sections extended infinitely far upstream or to the axis of the im- 
peller and thus yield no information concerning the flow behavior abead 
of and at the entrance to blades of finite thickness or blades which are 
not alined with the inlet stream. 

A rapid approximate method has been developed in reference 6 that 
will predict the blade surface velocities in centrifugal impellers. The 
accuracy of this method has been investigated for the region downstream 
of the inducer section. However, because of the lack of an exact solution 
of the flow in the inlet region, there has been no verification of the ap- 
proximate method there. 
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Consequently, a method for analyzing the flow in a centrifugal im- 
peller, including the inlet region, hy numerical solution of the partial 
differential equation governing the flow was developed at the NACA Lewis 
laboratory and is presented herein. This method was applied to a 48-inch- 
dlameter centrifugal impeller for both design and off-design flow con- 
ditions. The impeller is similar to that discussed in references 7 and 
8. The results of the exact solution are compared with the res\xlts of the 
rapid approximate method of reference 6. 


ANALYSIS 

The formulation of the problem and the proposed method of solution 
of the problem are discussed in this section. 


Statement of Problem 

The basic assumptions which are made concerning the physical nature 
of the flow determine the partial differential equation governing the 
flow. The assigning of proper boundary conditions to the problem then 
determines the particular solution of the p^tial differential equation. 

As8uBg)tions . - The flow is assumed to be steady, incompressible, and 
nonviscous. The assumption of steady, nonviscous flow is customary in 
compressor flow analyses. Several solutions have been obtained taking 
compressibility into account (e.g., refs. 1 and 2). The addition of the 
condition of compressibility complicates the solution by a factor which 
was considered to be out of proportion with its value in this case. 
Therefore, the fluid was assumed Incompressible. 

The further assumption is made that the flow is constrained to a 
blade -to -blade surface of revolution which is symmetrical about the im- 
peller axis. Although the flow is constrained to this surface, a varia- 
tion in the thickness of the stream sheet provides a closer approximation 
to the actual case. The shape of the stream surface in the axial -radial 
plane as well as the thickness variation is defined as a function of 
radial position which is specified at the beginning of the solution. 

The rear stagnation point is assumed to be located at the blade tip. 
The Kutta condition states that for a noncuspidate blade with a sharp 
trailing edge the rear stagnation point occurs at the tip. However, for 
an Impeller with a rounded trailing edge the location of the rear stagna- 
tion point cannot be predicted. It was necessary, therefore, to assume 
the location of the rear stagnation point. 

Differential equation . - In this analysis, the rl^t -handed cylin- 
drical coordinates r, 0, and z (see figs. 1 and 2) are used. All sym- 
bols are defined in the appendix. The angular velocity of the Impeller 
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Is denoted "by cd and the fluid density by p. The stream-sheet thick- 
ness in the z-direction is represented by b. The trace of the stream 
surface in the axial -radial plane is given by specifying z as a function 
of r. The slope dr/dz of this curve is equal to the tangent of the 
angle between the axis of rotation and the tangent to the trace of the 
stream surface in the axial-radial plane and is denoted by X (see fig. 

2). Thus, the resultant velocity w is given by 


w 


2 


w| + w| 



where the subscripts r and 6 indicate components in the r- and 
0-directions, respectively. 


The stream function 
equations : 


1 is defined by the following differential 


hi 


- bpWg 


(la) 


hi 

50 


rb pwp 


(lb) 


In this report all derivatives with respect to r shall be imderstood 
to mean derivatives with respect to r on the stream surface, that is, 
5/5r in this report shall correspond to the bold-faced h/hr of refer- 
ence 9, in which the differential equation for the type flow consid- 
ered herein is derived. With these definitions and assumptions, the 
differential equation of the flow becomes (see ref. 9, p. 35) 


hh 

5r2 


+ 



5 In b\ fn . _1_'\ 

~5F~;5? \ x^/ he^ 


2aSbp 


( 2 ) 


This equation, together with the boundary conditions, mathematically 
determines the problem. 


Boundary conditions . - This analysis of the flow is a boundary-value 
problem of the first kind or a Dlrichlet problem. Certain boundaries of 
the flow and the values of the stream function on these boundaries are 
specified. Furthermore, the flow is assumed to vary periodically in the 
circumferential direction, conq>leting a cycle in one pitch angle, the 
angular distance between two adjacent blade mean lines. The rotational 
speed of the Impeller and the weight flow through the compressor are also 
specified. The domain of the solution is extended sufficiently far up- 
stream and downstream, so that the flow is assumed to be uniform at the 
upstream and downstream boundaries. With the addition of these conditions, 
the problem is determined mathematically. 
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Method of Solution 


Superposition of four haslc solutions . - The dJJTferential equation 
(eq. (2) ) was' soiv-ed hy a superposition of four basic solutions. These 
four basic solutions form a set of linearly independent solutions such 
that all possible flows (including all tip speeds) are espre&sible as 
linear combinations of these basic solutions. The first of these, desig- 
nated tITq, is a solution of equation (2) with the condition that no flow 

crosses the upstream and downstream boundaries and cd » oDq ^ 0. 


The other three basic solutions, designated ijr^, and are 

solutions of the linear hamogeneous ecpiation obtained by equating the 
left side of equation (2) to zero. Thus, if 


L 



^ In b \ S i 
5r y 5r ^2 



£_ 

502 


then 1 ^^, -^fg, and liTg are solutions of 


L(ilf) = 0 


Because L is a linear operator, plus linear combinations of +1^ 
ijTg, and 1^2 will satisfy equation (2) for cs = 

m 

Boundary conditions for four basic solutions . - The flow region is 
represented by abcdef'gh in sketcia (a)'. 






e 



Sketch (a) 
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The upstream and downstreaifi houndaries, ah and de, respectively are 
placed sufficiently far from the "blades so that flow conditions can be 
assumed uniform at these stations. The angular distance from a to h and 
from d to e is one pitch angle. For all the basic solutions, the condi- 
tion that the flow is periodic about the axis of rotation with a period 
of one pitch angle makes it possible to obtain the solutions without a 
knowledge of the stream function along ab and gh. The finite-difference 
equation for points along these lines is obtained in the same manner as 
in reference 10. For the solution ilfQ in which the flow is that Induced 

only by the rotation of the Ingieller without any through flow, the value 
of along ah and de is specified zero, indicating no flow crossing the 
upstream and downstream boundaries. The values along the blade surfaces 
be and gf are also specified as zero. The solution to equation (2) for 
these boundary conditions is designated 

The through-flow solution, that is, the flow through the stationary 
blade row, is obtained from linear combinations of three basic solutions 
to the homogeneous equation obtained by equating the left side of equa- 
tion (2) to zero. All possible flows through the stationary blade row., 
can be represented by linear combinations of these basic solutions. It 
can be seen from the boundary conditions shown in the following table 
that these three basic solutions are linearly Independent. 


Basic 

solution 

Boundary values at 

Point a 

Point d 

Point e 

Point h 

■■ 

0 

0 

0 

0 


0 

0 

1 

1 


-1 

0 

1 

0 

■9 

-1 

-1 

0 

0 


The value of i{r along the trailing face is specified as zero and along 
the driving face as 1 for all three solutions. That these independent 
solutions are sufficient for the construction of nVl possible throu^ 
flows can be seen from the following consideration. 

Any flow is determined when the upstream and downstream flow condi- 
tions are specified. The velocity is assumed to be constant at stations 
ah and de. Hence, the stream function varies linearly from a to h and 
from d to e. Since the angular distance from a to h and from d to e is 
one pitch angle and the flow is assumed to vary periodically about the 
axis with a period of one pitch angle, the specification of conditions 
at points a, d, and h fixes the solution. The value of i|r is constant 
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along 'both the driving and trailing faces, with the difference between 
the function values being equal to The choice of 0 and ±1 

as the values of i|f at points a, h, d, and e and along be and gf repre- 
sents no restriction since these three basic solutions remain solutions 
to the homogeneous equation when changed by a multiplicative or additive 
constant . 

Coefficients of i|r^, ilr^, and in linear combinations. - The 

final solution f for any weight flow or rotational speed will be ob- 
tained from an equation of the form 


y = AqiItq + Aj^ti + ^2^2 + 

The coefficients Aq, A^_, A 2 } and A^ are determined by the specifica- 
tion of fo-ur Independent physical conditions: (l) the rotational speed, 

( 2 ) the weight flow, (3) the location of the rear stagnation point, and 
( 4 ) irrotatlonality of the absolute flow. 

The coefficient Aq is determined by the rotational speed and is 
given by 



( 4 ) 


That is, Aq is the ratio of the rotational speed co for the desired 
solution to that used in obtaining the basic solution cdq. 

The change in Y across one blade passage is equal to the weight 
flow through a single passage. Therefore, 

■^1 + ^2 + A 3 = M ( 5 ) 


where M is the weight flow through a single passage. 

The rear stagnation point is assumed to be at the blade tip. Thus, 
at the tip. 


s> 




Wg = 0 


or 



( 6 ) 
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This derivative is expressed in finite-difference form for the grid point 
at the blade .tip and 'with equation (3) yields a linear relation in Aq, 

Aq^^ Ag^ and '^ 3 * 


The absolute flov is Irrotational, so that^ if r^^ and rg are 
radial stations upstream of the blade row, the following equation holds: 


4 


tl 



^^ 9,2 + - 



(vg^l + asr^)r^de = 0 


( 7 ) 


where the subscripts 1 and 2 Indicate values along the lines r = r^ 
and r = rg, respectively. If r^ is chosen equal to the value of r 
at the upstream boundary, equation (7 ) becomes 



0,2 


d 6 = - 2 a>itr| 


( 8 ) 


since Wq ^ is equal to When the stream-function definition (eq. 

(la)) is introduced, equation ( 8 ) becomes 


^ d0 = 2 «upbgrg (9) 

Equation (9) can be Integrated nvunerically to yield 
A^, Ag, and A 3 . 

Equations (4), (5), (S), and (9) form a system 
Independent equations in four unknowns Aq, A^_, Ag, 

Mvimerlcal method of obtaining basic solutions , 
tion was covered with a network of grid lines whose 
grid or nodal points, as shown in figure 1 . 

The solution for a given set of boundary conditions of the differ- 
ential equation was obtained at each of these grid points by solving the 
set of linear simultaneous equations obtained when the differential equa- 
tion is written in finite-difference form for each grid point. A five- 
point system was used in the finite-difference approximation of the deriv- 
atives. This procedure is equivalent to approximating the stream function 
by a fourth-degree polynomial in the neighborhood of the grid point. The 


a linear relation in 


of four simultaneous 
and A 3 . 

- The region of solu- 
Intersectlons form 
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solution of the set of n linear simultaneous equations was obtained on 
high-speed digital computers hy the matrix method outlined in reference 
10. Since there were four basic solutions, four sets of n simultaneous 
linear equations were solved by this process. 


Numerical Example 

The previcfusly outlined method was applied in order to analyze the 
flow in a 48-inch-dlameter radial-inlet centrifugal impeller. A de- 
scription of the geometry of the impeller and the operating conditions 
for which the analysis was carried out follows. 

Geometry of impeller . - The impeller investigated was a 48-lnch-tlp- 
diameter radial-inlet centrifugal impeller having 18 blades, similar to 
that discussed in references 7 and 8. The sharp leading edge and blunt 
trailing edge were rounded as shown in figure 1 because of practical com- 
puting considerations. The blade coordinates are given in table I. The 
solution was obtained on the surface generated by rotating the mean blade- 
height line about the axis of rotation. This line was approximated by the i? 
following function; 


The streamline spacing in the axLal-radial plane is not known. Therefore, 
the stream-sheet thickness b in the z-direction was approximated by the 
blade height in the z-direction. This parameter was approximated by the 
following function: 


b = 0.07208 + 1.01517 (U) 

The parameter X is equal to dr/dz of the stream-surface trace in the 
axial-radial plane and from equation (lO) is given by 


1 ^ 0.041456 

^ (r - 0.40828)2 


( 12 ) 


Operating conditions . - Four solutions were obtained corresponding 
to four weight flows at a tip speed of 700 feet per second. These four 
weight flows, which correspond to those of reference 9, are shown in the 
following table; 


Case 

Weight flow, 
lb/ sec 

A 

14 

B 

26.25 

C 

32.10 

D 

44 
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RESULTS AM) DISCUSSION 

The results of the solutions obtained by the application of the pre- 
viously outlined method are presented in figures 3 and 4 , which show 
streamlines and constant relative velocity contours, respectively. Fig- 
ures 3 and 4 are projections on the r0 -plane, that is, the curvature of 
the stream surface in the axial -radial plane is neglected. The con 5 )ari- 
son with the results of the approximate method described in reference 6 
is made in figure 5. 

It was considered desirable to obtain the distribution of the veloc- 
ity at the blade nose in more detail than was obtained by this solution. 
Therefore, a refinement of the solution in this region was made by relax- 
ation methods and is presented in reference 11. 


Streamlines 

The distribution of stream function is shown by means of contours 
of constant stream-function ratio f/M in figure 3 for the four weight 
flows investigated. The Impeller tip speed was 700 feet per second for 
all four cases. 

Case A . - In figure 3(a) the streamlines for a weight flow of 14 
pounds per second are shown. This condition is the incipient surge 
weight flow for the experimental case (ref. 9). A large eddy attached 
to the driving face of the blade extends from r~1.30 to r~1.84 and 
almost one-third the distance across the passage between blades at its 
widest point. A stagnation point exists near the center of this eddy 
region. The major part of the flow is concentrated in the region near 
the trailing face, while the eddy and other relatively low-momentum fluid 
occupy half the channel. 

After the flow is turned to the blade direction in the inlet region, 
the flow outside the eddy remains roughly parallel to the trailing face 
of the blade to r ~1.5. Slip then begins to occur and the slip factor, 
the ratio of the mass-averaged absolute tangential velocity of the fluid 
at the tip to the absolute tip speed, is equal to 0.874. 

Case B . - The streamline pattern for the experimental peak pressure 
ratio and peak-efficiency weight flow (26.25 Ib/sec) is shown in figure 
3(b) . This weight flow is sufficiently high to eliminate the eddy on the 
driving face of the blade. However, a fairly large concentration of low- 
momentum air is still present, so that halfway through the impeller 50 
percent of the fluid occupies more than two-thirds the available flow 
area. Slip begins at r~1.50, and the slip factor is 0.873. 
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Case C . - The streamline pattern for a weight flow of 32.10 pounds 
per second (fig. 3(c)) Is similar to that for 26.25 pounds per second 
"because of the con5)aratlvely small change in weight flow. The flow con- 
tinues to shift toward the driving face, and the slip factor is 0.871. 

Case D. - In figure 3(d), the streamline pattern is shown for a 
weight flow of 44 pounds per second. In the investigation reported in 
reference 7 this represented the maximum weight flow. The flow is dis- 
tributed across the passage more nearly \iniformly than in the other ex- 
amples. The flow ceased to he perfectly guided at r~1.50, as occurred 
for all other weight flows. The slip factoiTfor this condition is 0.859, 
so that the total variation in slip factor from stxrge to maximum weight 
flow is only 0.015. 


Eelatlve Velocity 

As stated previously, the solution in the immediate nel^horhood of 
the nose was refined and is reported in reference 11. Hence, no discus- 
sion of the nose region is presented herein. 

Case A . - Contours of constant relative velocity ratio W (relative 
velocity divided, hy tip speed) are plotted in figure 4(a) for a weight 
flow of 14 po\inds per second. Downstream of the leading edge along the 
trailing face, there is a deceleration up to r <>'1.2. The velocity then 
remains constant to r~1.70, where a deceleration begins and continues 
to the stagnation point at the tip. Along the driving face, the flow 
decelerates after the slight acceleration just downstream of the stagna- 
tion point. This deceleration continues to the sta^snation point at 
r<-1.30, at which point the eddy begins to form. The eddy exists along 
the greater part of the blade. There is a slight acceleration near the 
tip. 


Case B . - In figure 4(b) the velocity contours for a weight flow of 
26.25 pounds per second are shown. After the slight deceleration just 
downstream of the leading edge, the velocity along the trailing face is 
constant (except for a small acceleration and deceleration at r ~ 1.3) to 
r ~ 1.70. The flow then decelerates to the stagnation point at the tip. 
Along the driving face, the velocity decreases to r ~ 1.36 and then is 
constant to r ~ 1.80. The flow begins to accelerate and then decelerates 
to the tip stagnation point. 

Case C . - Figure 4(c) presents the velocity contours for a weight 
flow of 32.10 pounds per second. The velocity along the trailing face is 
constant from r**1.0 to r~1.3 after the deceleration just downstream 
of the leading edge. At r- 1.3, a slight acceleration followed by a 
rapid deceleration occurs. This velocity peak is caused by the beginning 
of more rapid blade curvature at that point. This effect of the blade 
curvature is more evident for the higher wel^t flows than for the lower 
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iDecgoise of the absence of the deceleration on the trailing face just down- 
stream of the leading edge for the higher weight flows. Downstream of 
this rapid deceleration caused by the blade ciirvature, the velocity is 
constant from r^-l.SS to r~1.7. The usual deceleration to the stag- 
nation point begins at that point. On the driving face, the flow decel- 
erates from Just downstream of the leading edge to r ~1.4. The velocity 
is then constant to r~1.7 where the flow begins to accelerate slightly, 
followed by the deceleration to the stagnation point. 

Case D . - The velocity contours for a weight flow of 44 pounds per 
second are plotted in figure 4(d). The flow along the driving and trail- 
ing faces is similar to that for 32.10 pounds per second, because in both 
cases the inlet flow is directed toward the trailing face. The local ac- 
celeration and deceleration on the trailing face at r"1.3 is more 
pronounced. 

Angle of attack . - The approximate mean angle of attack, that is, 
the angle between the mean flow direction at the inlet and the tangent to 
the blade mean line, was computed from the rotational speed and the aver- 
age inlet velocity. The average inlet velocity was computed from the 
weight flow and the annular area. Two values were used for the annular 
area; (l) the total annular area with no blade blockage assumed, and (2) 
the total annular area minus the blockage caused by the blades. The 
thickness of the blades used in the latter computation was that at the 
1.4-foot radius, which was approximately the radius at which maximum blade 
thickness in the tangential direction occurred. The mean angle of attack 
across the passage at the 1.4-foot radius was also computed from the exact 
solution. These are compared in the following table: 


Weight flow, 
lb/ sec 

Angle of 

attack, based on - 

Unblocked 

annulus, 

deg 

Blocked 

annulus, 

deg 

Exact 

solution, 

deg 

14 

12.3 

9.2 

0 

26.25 

.1 

-4.6 

-5.3 

32.10 

-4.9 

-10.0 

-9.0 

44 

-13.6 

-19.0 

-15.4 


The sign convention for the angle of attack is such that a, positive angle 
of attack indicates flow directed towards the driving face of the blade. 
From the comparison of these angles of attack, it is apparent that the 
mean angle of attack is best predicted by basing the calculations on the 
annular area with blade blockage considered. The poor agreement between 
the mean angle of attack of the exact solution and that based on the 
blocked-inlet annular area at the lowest weight flow is probably caused 
by the ed(3y. It also appears from these average angles of attack that 
the original blade shape shown by the dashed line in figure 1 would per- 
form better at a weight flow of 26.25 poxmds per second. The blade angle 
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of the driving face Just downstream of the rounded leading edge Is 57°, 
whereas the angle between the mean line and the radial direction is 62°. 
Thus, at the weight flow of 26.25 pounds per second the average inlet 
flow angle would approximately equal the driving-face blade angle. 


Con^arison with Approximate Method 

The blade surface velocities for the four weight flows were computed 
by the rapid approximate method of reference 6 in order to determine the 
accuracy of the approximate method, especially in the inlet region, by 
conparison with the matrix solution; In figure 5, the surface velocity 
ratios as computed by the approximate method are compared with the sur- 
face velocity ratios as determined by the matrix method. The velocity 
ratios are plotted against distance s along the blade mean line. 

In the approximate method the mass -averaged velocity ratio Wjjj is 
assumed equal to the average of the blade surface velocity ratios 
The mass-averaged flow angle is assumed equal to the average of the 

blade surface angles Pg^^ for values of s<1.34. For s>1.34, Pj^ is 

approximated by a parabolic variation (with r) between the average of 
the blade surface angles at s = 1.34 and the flow angle at the tip as 
determined by the slip factor. The slope dPg^/dr at s = 1.34 is used 

as tbe third condition to determine the parabola. The approximate flow 
angle as given by the parabolic variation is also denoted by Pg^. In 

order to evaluate the validity of these assumptions, p^ is compared with 
pgy. and Wjj^ is compared with Wg^ in figures 6 and 7, respectively. 

It was inposslble to obtain meaningful results in the region of the 
roTinded leading edge. This failiire of the approximate method is caused 
by the invalidity of the assumption that Pj^ is equal to Pav* ^°^ 

values of s such that 0.07<s<0.6, the surface velocities are pre- 
dicted adequately for cases B and C but poorly for cases A and D. For 
case A the failure is probably caused by the change in the relative size 
of Pg^ and Pjj^. At s = 0.6, Pg^. is greater than Pjj but then becomes 

less than p^ as s decreases. This trend is not followed in the other 
cases. For case D the agreement of ¥av "with Wm and of Pav with Pm 
is poorer for s<0.6 than for the other cases. The trends for case D, 
however, were the same as for cases B and C. 

In the region 0.6^ s <1.8, the agreement between the approximate 
solution and the exact solution is adequate (except in the neighborhood 
of s = 0.8 for cases B, C, and D) for approximately predicting both 
the velocity and the velocity gradients for all four weight flows. Such 
good agreement is surprising for case A because the. agreement between 
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WJ^ and and "between and pg^^ is not nearly so good for case 

A as it is for the other cases; the eddy attached to the driving face for 
case A is the cause of this poor agreement. 

For s>1.8, the agreement between the approximate method and the 
exact method is poor and Pg^^ is less than Pj^ for a1 1 four cases ex- 
cept near the tip in case A. This reverse in trend, together with the 
fact that the "blade surface angles differ "by a large amount, pro"ba"bly is 
the reason for the failure of the approximate method near the tip. In 
the region s>1.34, pg^y is determined "by a parabolic approximation to 

Pjjj. Since Pg^ and Pjjj agree fairly well at s = 1.34 and at the tip, 

the failure of the paraholic approximation must "be attributed to in- 
correct values of d^g^/dr at s t= 1.34 or to the Inadequacy of a, para- 
bolic approximation. 


SUMMAEI OF RESULTS 

A method for the solution of the inconpressible nonviscous flow 
through a centrifugal impeller (including the inlet region) is presented 
and is applied to a 48-lnch-dlameter centrifugal impeller. Solutions 
were obtained for four weight flows ranging from incipient surge to max- 
imum as determined by actual in^eller tests. The blade surface veloci- 
ties obtained by these solutions were compared with those obtained by a 
rapid approximate method of analysis. The following results were noted; 

1. A large eddy formed on the driving face of the blade at the in- 
cipient surge weight flow "but was not present for the three higher weight 
flows . 


2. The slip factor varied from 0.874 to 0.859 as the weight flow 
Increased. 

3. For weight flows of 26.25, 32.10, and 44 pounds per second, a 
local acceleration followed by a rapid deceleration occurred on the 
trailing face of the blades at a radius of about 1.5 feet, that is, where 
the blade begins to curve more rapidly. 

4. The rapid approximate method waa adequate in the inlet region 
just downstream of the leading edge for weight flows of 26.25 and 32.10 
pounds per second. These weight flows corresponded to mean angles of 
attack of -5.3° and -9.0°, respectively. 
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5. The mean angle of attack was heat predicted "by basing the approx- 
imate computation on the weight flow, the tip speed, and the ann u lar area 
minus the blockage of the blades. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, March 4, 1955 
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APEEMDIX - SYMBOLS 


The follcwing syiribols are \ised in this report: 


coefficients in eq. (4) 


stream-sheet thickness in z-direction, ft 
weight flow through single passage, Ib/sec 
radial distance, ft 
distance along hlade mean line, ft 
ratio of relative velocity to tip speed 
relative velocity, ft/ sec 
axial distance, ft 

angular coordinate in relative system, radians 
slope of trace of stream sxxrface in axial-radial plane 
fluid density, Ib/cu ft 
stream function, eq. (l) 


basic solutions 


angular velocity of impeller, radlan/sec 


approximation to mass-averaged value used in approximate method 
mass -averaged value 
component in radial direction 
ing)eller tip 

component in tangential direction 

value in basic solution of eq. (2) 

conditions along r = r]_ and r = r 2 , respectively 


1,2 
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TABLE I. - MODIFIED BLADE COORDINATES 


Driving face 

Trailing face 

r 

0 

r 

0 

1.0405 

0.73321 

1.0285 

0.72685 

1.0521 

.72685 

1.0352 

.69779 

1.0696 

.70353 

1.0405 

. 68973 

1.0740 

.69780 

1.0518 

.66871 

1.1190 

. 63965 

1.0696 

. 64304 

1.1278 

. 62944 

1.0971 

. 61056 

1.1711 

.58150 

1.1278 

.57638 

1.2324 

. 52335 

1.1511 

. 55241 

1.2441 

. 51335 

1.2144 

. 49426 

1.3095 

.46520 

1.2441 

.47164 

1.3604 

.43582 

1.2974 

.43611 

1.4262 

.40705 

1.3525 

.40704 

1.4767 

.39101 

1.3604 

.40350 

1.5930 

.36851 

1.4283 

.37798 

1.7093 

.36080 

1.4767 

.36464 

1.8256 

.35976 

1.5600 

.34890 

1.9419 

.35976 

1.5930 

.34536 

2.0000 

.35497 

1.7093 

.33994 

2.0123 

.34890 

1.8256 

. 34161 



1.9419 

.34265 



2.0000 

.34432 



2.0123 

.34890 




18 


WA.CA TN 344:8 



Tlgure 1. - Grid STatea for oror-all aolutlon. 
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Plgfure 2. - Axial -radial plane view of 48-inch- 
diameter centrifugal impeller. 


X 

Axis of rotation 


second (case A) 


(a) Weight flow, 14 pounds per 

Figure 5. - Streamlines. 
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(Ti) Welglit flow, 26.25 pounds per second (case B) 
Figure 3. - Continued. Streamlines. 
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(d) Weight flov, 44 pounds per second (case D) . 
Figure 3. - Concluded. Streamlines. 
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(b) WelsSit flow, 26.25 potmde per second (case B) . 

Figure 4. - Continued- Contours of constant relative velocity ratio Vf. 
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(o) Weight flow, 32.10 pounds per second (case C) . 



Figure 4. - Continued. Contours of constant relative velocity ratio W 
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(d) Welgiit flow, 44 poundB per eecond (case D) . 

Figure 4. - Concluded. Contours of constant relative velocity ratio W. 
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1 

1 



Figure 5. - Continued. Congjarlaon of tlade-eurfaoe rcloeltf ratios obtained ty approximate and 
matrix methods. 
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(c) Weight flov, 32.10 pouzide par second (case C)» 


Figure 5. - Continued. Coiqparlson of hXade-surface 7eloclty ratios obtained by approximate and 
matrix methoda. 


c 

I 


i 


Blade-ourfaoe velocity ratio 
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(d) Weight flow, pounds per second (case D) . 

Figure 5. - Concluded. Comparison of blade-surface velocity ratios obtained by 
approximate and matrix methods. 
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(a) Weight flov, 14 pouMe per eecond (case A) . 


Figure 6. - CoD5>arlson of mass -aver aged flow angle of exact solution p_ 
and approximation of used In proximate method p^^^. 


r 

c 

c 

r 
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0 .4 .8 1.2 1.6 2.0 2.4 

Distance along blade mean line^ a, ft 


(c) Weight flow, 32.10 pounds per second (case C) . 

Figure 6. - Continued. Con^arlson of mass-averaged flow angle of exact 
solution and approximation of p^^ used In approximate method Pav 





t 


Distance along tlade mean line, m , ft 

(a) Weight flov, 14 pounds per second (case A) . 

Figure 7. - Con 5 >axison of mass-averaged relative velocity ratio W, 
average of hlade-surface velocity ratios Wgy. 
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Figure 7 
ratio 


.6 

Distance along 1310(16 mean line, s, 


ft 


(t) Weight flow, 26*25 pounds per second (case B) . 


► - Continued. Comparison of mass-averaged relative velocity 
Wm with average of hlade-«urface velocity ratios Waw* 




2.0 


2 . 


Figure 

ratio 


.6 

Distance along "blade mean line, a, ft 

(c) Weight flow, 32.10 pounds per second (case C) . 

. - Continued. ConjiaxlBon of maaB-avecaged relative velocity 
Wjji with average of hlade-Burface velocity ratios W . 




Dlstenee along blade mean line, s, ft 
(d) Wel^t flov, pounds per second (case D) . 


Figure 7. 
ratio ’ 


- Concluded. Con 5 >arison of mass-averaged relative velocity 
L -with average of blade~surface velocity ratios WoTr. 





